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The phase diagram of two-color QCD with non-zero chiral chemical potential is studied by means 
of lattice simulation. We focus on the influence of a chiral chemical potential on the confine¬ 
ment/deconfinement phase transition and the breaking/restoration of chiral symmetry. The simu¬ 
lation is carried out with dynamical staggered fermions without rooting. The dependences of the 
Polyakov loop, the chiral condensate and the corresponding susceptibilities on the chiral chemical 
potential and the temperature are presented. The critical temperature is observed to increase with 
increasing chiral chemical potential. 
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I. INTRODUCTION 

For the vacuum state of QCD as well as for the properties of finite-temperature QCD the existence of non¬ 
trivial topological excitations is important. Well known are instantons [T] as classical solutions in Euclidean 
space as tunneling events between topologically different zero-temperature vacua. The role of topology for the 
solution of the famous Ua{^) problem has been recognized very early [3[S]. Presently, the anomalous breaking 
of the C/yi(l) symmetry above the deconfinement transition/crossover is under intense investigation in the lattice 
community (see, e.g.. Ref. a)- 

It is known by now from lattice QCD at zero temperature that a (fractal) low-dimensional (laminar) topo¬ 
logical vacuum structure is discernible at very fine resolution scale [S], while localized instanton-like structures, 
actually prevailing at an infrared scale, are believed to explain chiral symmetry breaking mm- 

The gluon fields contributing to the path integral at finite temperature correspondingly may contain 
calorons mm- Adapted to the non-trivial holonomy they have a richer structure (in terms of “dyonic” con¬ 
stituents) than instantons. The changes of this structure at the QCD phase transition are presently under 
study [TOl ITT] . 

Some time ago the gluonic topological structure and the famous axial anomaly have been proposed to be 
immediately observable (and controllable) through the generation of P and CP violating domains (violating 
also translational invariance) in heavy ion collisions [mill]. It has been demonstrated by detailed numerical 
calculations [HHll that macroscopic domains of (anti)parallel color-electric and color-magnetic field can emerge 
in a heavy ion collision creating an increasing chiral imbalance among the quarks which are deconfined due to 
the high energy density. In this situation, the magnetic field created by the spectator nucleons may initiate a 
charge separation relative to the reaction plane (parallel to the electro-magnetic field) [H] . The resulting charge 
asymmetry of quarks would become observable in terms of recombined hadrons (chiral-magnetic effect) [161117) . 
The strength (and particularly the dependence on the collision energy) of this effect has been theoretically 
studied and proposed to be a signal of the transient existence of liberated quarks [Hiiaiii]- 

In recent years the dependence of the chiral and deconfinement transitions on the magnetic held has been 
investigated both in models and ab-initio lattice simulations, see e.g. [niin]. It remains an open question 
whether the phase transition from quarks to hadrons, i.e. the onset of conhnement and chiral symmetry breaking 
(and vice versa), depends on the chiral imbalance. 

In this paper we study the change of the phase structure induced by a chiral chemical potential in an 
equilibrium lattice simulation. We mimic the topological content (of a CP violating, topologically nontrivial 
gluonic background in a heavy ion collision event) by inducing a chiral imbalance, which is provided (i.e. frozen) 
by a non-zero chiral chemical potential. In this setting, the modihcation of the phase diagram by the chiral 
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chemical potential /i 5 has been studied mainly by means of effective models [HHiT], the predictions of which 
will later be compared with our results. 

On the lattice, contrary to the case of non-zero baryon chemical potential, simulations with non-vanishing /Z 5 
are not hampered by a sign problem^. Thus, one can employ the standard hybrid Monte Carlo algorithms. Such 
lattice simulations with fj .5 ^ 0 were already performed in Ref. [551 [53]. The main goal of these papers, however, 
was the study of the chiral magnetic effect. Therefore, the phase diagram was not systematically studied. 

In Refs. [SOj [3T| we have carried out the first lattice study of the phase diagram with non-zero chiral chemical 
potential. It was performed in SU{2) QCD with four flavors, which we have considered as a simplified model 
of QCD. In this article we extend this investigation generating considerably more data and performing a more 
detailed analysis. In this paper we do not attempt to analyse the topological structure that would reflect the 
presence of the chiral chemical potential. 

One reason for choosing the SU{2) gauge group is that less computational resources are required for this pilot 
study than for full QCD. The second reason is that we have already carried out two-colour QCD computations 
with an external magnetic field |32[ 13,1j . Furthermore, the four flavor case results from our choice of staggered 
fermions as lattice regularization while we avoid to take the root of the fermion determinant, which would allow 
to reduce the number of flavors. The “rooting procedure” might introduce further systematic errors at finite 
lattice spacing. 

In section II we introduce the model and its lattice implementation and define the quantities we measure. 
In section III we present our results, and section IV is devoted to their discussion and to the formulation 
of conclusions. In the Appendix we discuss the question of renormalizations refering to explicit analytical 
calculations in perturbation theory. 


II. DETAILS OF THE SIMULATIONS 


We have performed simulations with the SU{2) gauge group. We employ the standard Wilson plaquette 
action 

( 1 -]^Trt/^^(x)y ( 1 ) 

For the fermionic part of the action we use staggered fermions 


Sf = 

X 

+ ^ X! - '4>xUl{x)tljx+^)+ (2) 

Xfl 

T ^ s{x')(^'lpx+6Ux-\-S^x'4^x '^xU: 

X 

where the are the standard staggered phase factors: rii{x) = 1 , 77 ^( 0 ;) = ^ _ 2 , 3 , 4 . 

The lattice spacing is denoted by a, the bare fermion mass by m, and 7^5 is the value of the chiral chemical 
potential. In the chirality breaking term s{x) = (—1)“^, <5 = (1,1,1,0) represents a shift to the diagonally 
opposite site in a spatial 2^ elementary cube. The combination of three links connecting sites x and x + 5, 


Ux-\-5,x 


Ui{x + Cj + ek)Uj{x + ek)Ukix) 

2,j,/c=perm(l,2,3) 


( 3 ) 


is symmetrized over the 6 shortest paths between these sites. In the partition function, after formally integrating 
over the fermions, one obtains the corresponding determinant. As mentioned above, we do not take the fourth 
root of this determinant in order to represent each flavor (“taste”) independently of the others. Thus, the 
continuum limit of our model corresponds to a theory of four (degenerate) flavors. 


^ In SU (2) there is no sign problem even in presence of a baryon chemical potential. 
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In the continuum limit Eq. ([^ can be rewritten in the Dirac spinor-flavor basis [331133 as follows 

r 

Sf = J d'^x^ + m + ^ 57574 )%- (4) 

i—1 

We would like to emphasize that the chiral chemical potential, introduced in Eq. corresponds to the 
taste-singlet operator 7574 0 1 in the continuum limit. 

It should be also noted here that the usual baryonic chemical potential after the discussion in Ref. [36], 
and also the chiral chemical potential as used in Ref. [29]) are introduced to the action as a modification of 
the temporal links by corresponding exponential factors in order to eliminate chemical-potential dependent 
quadratic divergencies. For staggered fermions this modification can be performed as well for the baryonic 
chemical potential. However, for the chiral chemical potential such a modification would lead to a highly non¬ 
local action [33] ■ Therefore, we decided to introduce ^5 in Eq. (H in an additive way similar to the mass 
term. 

It is known that the additive way of introducing the chemical potential leads to additional divergencies in 
observables. In the Appendix we present analytical and numerical investigations of additional divergencies in 
the Polyakov loop and the chiral condensate. Our study shows that there is no additional divergency in the 
Polyakov loop, whereas there is an additional logarithmic divergency in the chiral condensate. The latter is 
numerically small and does not effect the results of this paper. 

We have performed simulations with two lattice sizes N-r x = 6 x 20^, 10 x 28^. The measured observables 
are 


• the Polyakov loop 


L = 


1 


Nr 

E (Trll 174(ni, 712,743, 714 )), 

ni,n 2 ,n 3 714 — ! 


( 5 ) 


• the chiral condensate 




I 1 d 


NrN^ 4 d{ma) 


\ogZ = 


1 I 


W7V3 4 
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D 




• the Polyakov loop susceptibility 


( 6 ) 


Xi = ((^2) _ (i)2) ^ 


( 7 ) 


• the disconnected part of the chiral susceptibility 


\disc 


1 1 


(((Tr 


D 


-f) - (Tr 


I 


D -f 777,0 


n- 


( 8 ) 


The Polyakov loop and the corresponding susceptibility are sensitive to the confinement/deconfinement phase 
transition, whereas the chiral condensate in principle responds to chiral symmetry breaking/restoration. 

The simulations have been carried out with a CUBA code in order to perform the simulations using the 
Hybrid Monte Carlo algorithm on GPU’s. 

The parameters of our lattice calculation are the inverse of the bare coupling constant /3, the bare mass ma 
and the bare chiral chemical potential a/is (both in dimensionless units) as well as the lattice size. The physical 
temperature and the volume are given by 

U = (iV.a(/3))3 

^ 1 ( 9 ) 

a{P)Nr' 

To perform the scale setting (calibration) of the lattice, we use the results of Ref. [32]. There, the static 
potential in the same theory has been measured at zero temperature for several values of /3. The Sommer 
parameter ro was calculated in lattice units and compared to its physical value, which was supposed to be the 
same as in QCD, i.e. ro = 0.468(4) fm. It was found that the scaling function a{(3) in the region were we 
perform our measurements is well described by the two loop /3-function. Thus, for given /3 we can obtain e.g. 
the temperature T in units of MeV. For more details, see Ref. [33] . 
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FIG. 1: Polyakov loop and chiral condensate versus T for five values of Lattice size is 6 x 20^, fermion mass is m « 12 
MeV. Errors are smaller than the data point symbols. The curves are to guide the eyes. 


III. RESULTS OF THE CALCULATION 

We first performed simulations on a lattice of size 6 x 20^ for five fixed values of = 0,150,300,475,950 
MeV and different values of T. The fermion mass was kept fixed in physical units m « 12 MeV(m 7 r ~ 330 
MeV). The expectation values of the Polyakov loop and the chiral condensate are shown in Fig. The sharp 
change of the observables as functions of T indicates the onset of the deconhnement and the chiral restoration 
phase transition. It is seen that the temperature of both phase transitions increases with the chiral chemical 
potential. One also sees that the phase trasition becomes sharper for increasing chiral chemical potential. 
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FIG. 2: Polyakov loop susceptibility and chiral susceptibility versus T for three values of ^5 = 0,475,950 MeV. Lattice 
size is 6 X 20®, fermion mass is m « 12 MeV. In order to avoid a complete superposition of data points belonging to 
different ps values we applied a tiny shift along the T axis. 

To study the change of the critical temperature more quantitatively, we also calculated the chiral and the 
Polyakov loop susceptibilities. In order to make the figures readable we plot separately the susceptibilities for 
values ^5 = 0,475,950 MeV in Fig. and for values (15 = 0,150, 300 MeV in Fig. We see that increasing 
the value of the chiral chemical potential moves the position of the peaks of the chiral and Polyakov loop 
susceptibilities to larger values of f3. This means that the transition temperature increases. Our results do not 
show any splitting between the chiral and the deconhnement transition. We have htted the data for the chiral 
susceptibility near the peak with a gaussian function x = oq + oi exp ^ and extracted the critical 

temperature Tc{/3c)- The resulting dependence of the critical temperature on the value of the chiral chemical 
potential is shown on the Fig. |^and Table |T] The values of j3c and Tc and their uncertainties are calculated 
from a £t of 5-6 points in the vicinity of the peak by the Gaussian function given above. It should be noted 
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FIG. 3: Polyakov loop susceptibility and chiral susceptibility versus /3 for three values of fis = 0,150, 300 MeV. Lattice 
size is 6 X 20^, fermion mass is m 12 MeV. In order to avoid a complete superposition of data points belonging to 
different fis values we applied a tiny shift along the T axis. 


that for small values of the dependence of Tc on fj,^ is well described by the function Tc = a + 6 ^ 5 , but the 
larger ^5 is the larger is the deviation from this simple formula. 
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FIG. 4: The dependence of the critical temperature on the value of the chiral chemical potential. Lattice size is 6 x 20^, 
fermion mass is m 12 MeV. The curve is to guide the eyes. 


At5 (MeV) 

P. 

Tc (MeV) 

0 

1.7975 ± 0.0005 

195.8 ± 0.4 

150 

1.7984 ± 0.0009 

196.7 ± 0.7 

300 

1.8012 ± 0.0007 

198.9 ± 0.5 

475 

1.8116 ± 0.0003 

207.5 ± 0.2 

950 

1.8404 ± 0.0002 

233.3 ± 0.2 


TABLE I: The critical temperature Tc and the lattice parameter/3c as a function of the chiral chemical potential ^5 
obtained from the fit of the chiral susceptibility near the peak by a gaussian function x~ T a\ exp ^ ■ 

In addition to the calculations on the 6 x 20^ lattice we carried out simulations on a larger lattice of size 
10 X 28^. Note that this allows us to investigate larger values of /X 5 . The susceptibilities require large statistics, 
and therefore our current computational resources do not allow us to measure them on the larger lattice. In 
the simulations we kept the physical fermion mass fixed at m = 18.5 MeV {m^ ~ 550 MeV). 

On this lattice we also calculated the observables L and as functions of T for different values of 
The results for the Polyakov loop and the chiral condensate are presented in Fig. They underpin the fact 
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FIG. 5: Polyakov loop and chiral condensate versus T for four values of ps and lattice size 10 x 28^, fermion mass is 
ma « 18.5 MeV. Errors are smaller than the data point symbols. The curves are to guide the eyes. 


that the critical temperature grows when one increases The value of the chiral chemical potential used in 
our simulations was rather large (up to /X 5 = 3345 MeV), but we have not seen signals of a hrst-order phase 
transition, although the transition becomes sharper when the value of the chiral chemical potential grows. 

It is interesting to study how the observables L and (V'V') depend on /i 5 for fixed temperature. To do this 
we have measured the Polyakov loop and the chiral condensate as functions of ^,50 for three different values of 
P, P = 1.87 (158 MeV), 1.91 (186 MeV) and 1.95 (219 MeV), which for a vanishing chiral chemical potential 
corresponds to a temperature below the transition, in the transition region, and in the high temperature phase, 
respectively. The results of these measurements are shown in Fig. As can be seen from the figure, in the 
confinement phase the Polyakov loop remains almost constant with increasing chiral chemical potential. It 
means that if the system was in the confinement phase at ^5 = 0, it remains confined at /is > 0. Moreover, 
we observe the Polyakov loop to drop down with increasing /is both in the deconfinement phase and in the 
transition region. Thus, the system goes into the confinement phase for sufficiently large /ts. With other words, 
we conclude that the critical temperature increases with an increasing chiral chemical potential in agreement 
with our results obtained on the smaller lattice. It is worth mentioning that the behavior described above looks 
quite similar to the behavior obtained for two-color QCD in an external magnetic field [ini Eg. 
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FIG. 6: Polyakov loop and chiral condensate versus ps for three P values and lattice size 10 x 28^, fermion mass is 
ma « 18.5 MeV. Errors are smaller than the data point symbols. The curves are to guide the eyes. 


In order to study the dependence of the observables on the value of the quark mass for different values of 
chiral chemical potential, we chose four values of /3: one in the confinement phase, P = 1.8 (corresponding to a 
temperature T = 119 MeV), two in the transition region, P = 1.9 (T = 178 MeV) and P = 2.0 (T = 268 MeV), 
and one in the deconhnement phase /3 = 2.1 (T = 404 MeV). The lattice size used here is 10 x 28^. We chose 
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FIG. 7: Chiral condensate and Polyakov loop versus ma for three fis values at /I = 1.80 on a lattice of size 10 x 28®. 
P = 1.80 corresponds to the confinement phase. 


four values of the fermion mass ma and three values of the chemical potential /is = 0 MeV, /X 5 = 1115 MeV 
and /is = 2230 MeV. The expectation values of the Polyakov loop and the chiral condensate are presented as 
functions of ma for different /3 in Fig. 0 Fig. in Fig. 0 and Fig. [TO] 
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TABLE II: The parameters for the fits fi{ma) Eq. (101 and f 2 {ma) Eq. ( |11[ ) allowing to extrapolate to the chiral limit 
for P = 1.80 and /3 = 1.90 and various values of chiral chemical potential. The fit curves obtained with /i are shown in 
the left panel of Fig. 0111 


In Fig. 0we show the results in the confinement region at f3 = 1.80 (T = 119 MeV). The Polyakov loop is 
small and does not show any nontrivial behaviour. The chiral condensate remains almost constant when the 
value of the fermion mass changes. 

To extrapolate to the chiral limit (ma = 0) we performed chiral extrapolations in two ways. Since the value of 
the temperature is not very far from the transition temperature, we suppose that we can use the formula for the 
reduced three dimensional model [37] ■ See also [SShdO] . In that case we may use the ansatz >= /i(ma), 
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FIG. 8: Same as Fig. 7, but for /3 = 1.90 corresponding to the transition region. 
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FIG. 9: Same as Fig. 7, but for (3 = 2.00 corresponding to the transition region as well.. 
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FIG. 10: Same as Fig. 7, but for /3 = 2.10 corresponding to the deconfinement phase. 


where 


fi{ma) = ao + a\^ma + a 2 ma (10) 

with the non-analytic term coming from the Goldstone bosons. Such a parametrization has been used in m in 
the context of the finite temperature transition in QCD. As an alternative we also use the chiral extrapolation 
relevant for zero temperature, namely < ijjip >= / 2 (ma), where 


/ 2 (ma) = bo + bimalogma + b 2 ‘ma. 


( 11 ) 


The fits performed with the Eq. ( |Io| ) are shown as lines in Fig. ^ The corresponding parameters of both are 
given in Table [ll| As can be seen from table II and the corresponding figure, the chiral condensate extrapolates 
to a non-zero value in the chiral limit. 

At larger values /3 = 1.90 (T = 178 MeV), P = 2.00 (T = 268 MeV), and /3 = 2.10 (T = 404 MeV) and 
at zero chiral chemical potential the system is in the chirally restored phase, when m goes to zero. However, 
at /? = 1.90 (r = 178 MeV) and p-s = 1115,2230 MeV, the chiral condensate has non-zero expectation values. 
We have performed the same extrapolation technique as for = 1.80 and present also these results in Table [ill 
showing that for these values of the data are consistent with a non-zero values in the chiral limit. It implies 
that also in the chiral limit a non-zero chiral chemical potential shifts the position of the phase transition to 
larger temperatures. The plot for the Polyakov loop (right panel of Fig. confirms this observation: for the 
two smaller values of the fermion mass ma the values of Polyakov loop for ^,50 = 0 is several times larger than 
for non-zero values of so that the system for non-zero /is is deeper in the confinement region than for zero 
chiral chemical potential. 

At /3 = 2.00 (T = 268 MeV) and P = 2.10 (T = 404 MeV) the chiral condensate goes to zero in the chiral 
limit (Fig. lOl, thus the system is in the chirally restored phase. The dependence of the chiral condensate on 
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the value of the mass is almost a linear function, apart from the behavior for j3 = 2.00 (T = 268 MeV) at the 
largest chiral chemical potential under our investigation, fj,^ = 2230 MeV. This behaviour can be explained if 
one assumes that this point is near the phase transition. It is known that increasing the mass shifts the position 
of the transition to larger temperatures |32j . Thus, at smaller masses the system is in the chirally restored 
phase, while increasing the mass lets the system undergo the chiral phase transition. Again, the plot of the 
Polyakov loop confirms this suggestion. The Polyakov loop is small at larger masses implying that the system is 
in the confinement phase and increases when the mass is decreased, i.e. the system becomes deconfining. This 
behavior suggests that at the largest value of /i 5 = 2230 MeV the transition happens at even larger temperature 
than we have investigated. These results are in total agreement with the results above. 


IV. DISCUSSION AND CONCLUSION 

In this paper we have presented an investigation of the phase diagram of two-color QCD with a chiral 
chemical potential using lattice simulations with dynamical staggered fermions without rooting, i.e. four flavors 
in the continuum limit. We have calculated the chiral condensate, the Polyakov loop, the Polyakov loop 
susceptibility and the chiral susceptibility for different values of the temperature T and the chiral chemical 
potential on lattices of size 6 x 20^ and 10 x 28^. The main result is that at non-zero values of the chiral 
chemical potential the critical temperatures of the confinement/deconfinement phase transition and of the chiral- 
symmetry breaking/restoration phase transition still coincide and that the common transition shifts to larger 
temperatures as p-s increases. It is shown that this conclusion remains true in the limit of zero quark mass. 
This result is in contradiction to the results within different effective models of QCD piHnniEiiiEii where the 
critical temperature is said to decrease as /is increases. 

We concede that we study SU(2) QCD with Nf = 4 quarks which is different from what is considered in the 
effective models, where SU(3) QCD with Nf = 2 quarks is considered. It would therefore be very interesting 
for a closer comparison to examine the results of these models by a lattice simulation. 

It is likely, however, that the contradictions are rather a consequence of the fact that the critical temperature 
is not a universal parameter but crucially depends on the structure and parameters of the effective models. It 
is unclear to what extent they describe the actual behavior of finite temperature QCD. It should be also noted 
that some predictions obtained in different effective models contradict among each others. For instance, the 
dependence of the chiral condensate on /is obtained in papers inilll] is opposite to the result obtained in paper 
pS] . The latter paper claims that the chiral condensate grows with the chiral chemical potential, which is in 
agreement with our results. 

In addition to the dependence of the critical temperature on the chiral chemical potential, some effective 
models predict that - beginning from some critical value of /Xg (/Xg « 400 MeV in Ref. [H], /i§ = 50 MeV 
in Ref. |22] 1 - the transition from the hadronic to the plasma phase turns into a first order transition. In 
our simulations we see that the transition becomes sharper as we increase /xg, but we don’t see a first order 
phase transition up to /xg = 3345 MeV, which would manifest itself as discontinuity in Polyakov loop and chiral 
condensate at a sufficiently large spatial extension of the lattice. In this connection a finite size scaling analysis 
would be valuable but is outside the scope of this paper. 

Besides within effective models, the phase diagram of QCD in the (/xg, T)-plane was studied in papers 
in a framework of Dyson-Schwinger equations. The authors of these papers found that the critical temperature 
rises with /xg and the “phase transition” actually is always a crossover. These results are corroborating the 
results of our paper. 

We would like also to mention the paper |44j . In this paper the authors address the question of universality 
of phase diagrams in QCD and QCD-like theories through the large-iVc equivalence. Using the results of this 
paper one can show that at large-Vc and the chiral limit there is equivalence between QCD phase diagram at 
finite /xg and QCD phase diagram at finite isospin chemical potential /x/. The chiral condensate of QCD with 
/Xg ^ 0 can be mapped to the pion condensate of QCD with /x/ ^ 0. In latter theory the pion condensate and 
critical temperature T^. of pion condensation increase with /x/. Despite the fact that we considered N^. = 2 one 
can expect that chiral condensate and Tc in our theory also increase with with /xg. We believe this is one more 
fact in favour of our results. 
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Appendix A: Ultraviolet divergences in the chiral condensate 


The fermion propagator including the chiral chemical potential for ’’naive” lattice fermions can be written in 
the following form 




^a/3 


EE' 


^ip(x-y) 


{p} 


-* IpSiniPf,) +m + HBlolb 
sin'^{po) + (|p| — s/rs)^ -I- 


X P{s), 


(Al) 


Pis) = 


jisin{p,) 


\P\ 


-7o75 


* = 1,2,3, 


|p|2 = sin^{pi) + sin^{p2) + sin^ips), 

Pi = i = l,2,3, ni = 0,...,Ls-l, 

2tt tt 

P4 = + 714 = 0, ...,Lt - 1. 

Here m and ^5 are mass and chiral chemical potential in lattice units, a, j5 are color indices, the sum is taken 
over all possible values of ( 711 , 712 , 713 , 714 ) and s = ± 1 . 

In the limit —)■ 00 the condensate per one fermion flavour can be written as 


d'^p 


V- 

f ^ air 


1 


('0*/’) a;)] ^ ^ sin'^ipo) + (|p| - ’ 


(A2) 


To calculate the integral in formula (A2| we use the standard approach, which separates the main divergence 
from the rest of the integral 


d'^p 


V- 

Z —^ 07 


r-Ti- j4 


d^p 


(27r)4 41^ siiPipo) + (IpI - 5 / 45)2 -I- m? sin‘^{po) -I- \p\^ + m? 


(A3) 


fir j4 


d'^p 


E 


(27r)4 ^\sin'^{po) + (|p| — 5 / 75)2 -|- 7772 sin‘^{po) + Ipp m? J 


The first term in this expression is just the loop integral for (t/’V') without chiral chemical potential. The 
expression for this integral in the continuum limit a —>■ 0 can be found in [45] 


d'^p 


(27r)4 57772 (po) -I- |p|2 + m? 


= 0.619734 -f 777 " 


0.345071)+ 0 (m->). 


V 




(A4) 


The second term in (A3) is proportional to /tg and contains a logarithmic divergence in the continuum limit. 
This divergence can be calculated using saddle-point method. To calculate the second term in ( |A3[ ) we subtract 
the leading divergence and calculate the rest of the integral in the limit a —^ 0. The result of the calculation 
can be written as follows 


r d^P ^ 

(2^)4 


1 


1 


sirP{po) + (IpI — 5 / 75)2 + m? sirP{po) + jpp -|- 7772 

log{m ) ^ Q 57 j^Q 35 ^ _|_ 0 {m*, pl\). 


= -4 


(A5) 


We have checked the last formula numerically. Introducing a physical mass vripa = m and a physical chiral 
chemical potential one can write the expression for the condensate in physical units as 


(7/,7/7)p = 0.309867 


777„ 


; f log{mpa) 


2 + ^47 I 

72 y TT^ 


^ - 0.172536^ -f 777p(/7f)2 + 0.167759^ . (A6) 
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From the last formula one sees that at one-loop level the inclusion of a non-zero chiral chemical potential leads 
to an additional logarithmic divergence. We believe that this conclusion persists if one takes into account 
radiative corrections to formula (A6). To see this, consider some Feynman graph of radiative corrections to 
formula (A6). The expansion of this graph in powers of ^5 can be reduced to the inclusion of a dimension-3 
operator to the graph that diminishes the power of divergence by one unit per one power of fj,^. Expansion of 
the chiral condensate in powers of fj,^ contains only even powers. The main divergence of the chiral condensate 
is ~ 1/0^. So, the next to leading term in /is expansion is two powers smaller than 
logarithmically divergent. 


1/a^, i.e. it is at most 


lattice results 
divergence 


0.07 


0.08 0.09 


0.1 


0.11 


0.12 


0.13 


aT 


FIG. 11: The difference between chiral condensates at /is = 0 MeV and /is = 330 MeV for different lattice spacings 
(T = 270 MeV). Square points correspond to lattice measurements, circle points correspond to the contribution of 
logarithmic divergence (A 61 to the difference. 


Further it is important to understand how the additional divergence connected to non-zero /is effects the 
results of this paper. To estimate this we fixed the physical values of temperature T = 270 MeV, quark mass 
nip = 33 MeV and measured the condensate in the deconfinement phase for /ts = 0 and /is = 330 MeV for 
different values of lattice spacing a. In particular, we took the following lattice parameters: 16^ x 8 /3 = 1.9500, 
20^ X 10 /3 = 2.0047, 24^ x 12 /3 = 2.0493 and 28^ x 14 /3 = 2.0870. The difference between the measured 
chiral condensates at zero and non-zero /is for different lattice spacings is shown in F ig|ll[ In addition to 
the lattice measurement, in Figwe plot the contribution of logarithmic divergence (A6) to the difference 

From Fig[^ one sees that uncertainty of the calculation doesn’t allow to confirm the logarithmic type of 
the divergence. However, the variation of the lattice results with variation of the lattice spacing is very slow 
what allows us to assume that there are no other divergences different from the logarithmic one. Note also that 
the value of the /is 7^ 0 contribution to the chiral condensate obtained from the one-loop expression for the 
logarithmic divergence is by a factor 2-3 smaller than the lattice results. The difference can be attributed to 
a non-zero temperature effect and radiative corrections. This fact allows one to state that the renormalization 
effect of the one-loop expression for the logarithmic divergence (A6) is not very large. So one can use it to 


estimate the contribution of ultraviolet logarithmic divergence in the confinement phase. 

The calculation shows that the characteristic values of the difference (1^1 /') — ('i/'' 0 )/j 5 =o in the confinement 
phase for different values of lattice parameters used in this paper are approximately by one order of magnitude 
larger than the additional ultraviolet logarithmic divergence due to /is 7^ 0 estimated according to formula 
(A6). For instance, from Fig. f^one can see that at temperature T = 200 MeV (close to the phase transition 
at fj ,5 = 0) the difference (('^Vv^g^aooMeV ~ fj. 5 =o )~ 2. At the same time the additional ultraviolet 

divergence according to formula (A6I gives {{'4’tp)add)/T^ ~ 0.1. Note also that the position of the phase 
transition manifests itself as a peak of the susceptibility. Evidently, there is no peak due to the additional 
ultraviolet divergence. All this allows us to state that the conclusions obtained in this paper are not affected 
by the additional ultraviolet /ts 7^ 0 divergence in the chiral condensate. 
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Appendix B: Ultraviolet divergences in the Polyakov loop 


At the leading order approximation in the strong coupling constant g the Polyakov loop can be written in 
the following form 


L = l-P^- 


'drU-(r,d) = l-4 / ^ 




(Bl) 


where DQQ{t,x) and DQQ{qo,q) represent the propagator of the temporal components of the gluon field in 
coordinate and momentum space, correspondingly, f3 is inverse temperature, and the summation over color 
index a is assumed. From eq. (Bl) it is clear that additional divergences connected to non-zero at one-loop 


level can appear from the fermion self-energy part of the gluon propagator. The study of the divergences with 
’’naive” fermions in the self-energy part of the propagator is cumbersome. For this reason we use momentum 
cut of regularization procedure to study the divergences in the Polyakov loop. 



FIG. 12: Renormalized Polyakov loop (B4| as a function of lattice spacing (T = 270 MeV). 


The fermion contribution to the one-loop self-energy of the gluon propagator D^q can be written in the 
following form 


noo(^) — 


4 ^ 

Sl,S2 ' 


d^k ^ _ -fcg -I- (^5 - Si|fc|)(/r5 - S2|fc + ^) + 

(27!')'* (/c2 _|_ (ifci _ + m?){kl + {\k + ^ - S2gbY + 


1 -|- S1S2 


k- [k + q) 

\k\\k + ^ 


(B2) 


where q is the external momentu m, si , S2 = ±1. 

To proceed we expand equation (B2) in ^5. Evidently each new term in this expansion diminishes the power of 
divergence by two units. That means that beginning from /r| term the series is ultraviolet convergent. The first 
term (noo(^|/25=o) is just the expression for the self-energy without chiral chemical potential. So, to calculate 
the divergence which result s fro m nonzero /is one needs to study only the divergence of the second term. Taking 
the derivative of equation (B2) with respect to we get 


noo(^ - noo((7)|/i5=o — Vs J 


d^k 


1 


(27r)‘^ \ {k^ + |fcp -I- ■m?){kQ + \k + + vt?) 


(B3) 


2m^ -I- 4|fc + q\^ + &k ■ {k -\- (tj — 2fcg 


8|fcp(TO^ -\-k ■ {k + q) — fcg) 


(fcg -I- |A:p -I- m?Y{kQ + \k + (f\^ + m?) (fcg -I- |A:p -I- w?Y{kQ + \k + (f\^ + w?) 

4(m^ — fcg)(fc ■ {k + q)) + 4\k\'^\k + 


(fcg -I- |fcp -I- nn?Y{kQ + \k + q\^ + im?Y 


Oigt) 


One can easily check t hat this expression is free from ultraviolet divergence. Moreover, a calculation shows that 
at large (f expression (B3) behaves as noo(^ —noo(g)|,i5=o ~ log-I-C 2 ), where ci,C 2 are some constants. 
So, at one-loop level there is no ultraviolet divergence in the Polyakov loop which results from non-zero chiral 
chemical potential. Similarly to the previous section, one can argue that there is no ultraviolet divergence in 
the g 5 ^ 0 contribution to the Polyakov loop at higher loops. 
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In order to numerically study the role of divergences in the Polyakov loop due to non-zero chiral chemical 
potential, we used the following definition of renormalized Polyakov loop 

^ren ~ y : 

where L is the Polyakov loop and is the Polyakov loop at the same lattice but with /X5 = 0. Evidently 

Lre„ doesn’t contain divergences that are usual for simulations with zero chiral chemical potential. Similarly 
to the previous section we fixed the physical values of temperature T = 270 MeV, quark mass nip = 33 MeV 
and measured the L^en in the deconfinement phase for ^,5 = 330 MeV and different values of lattice spacing a. 
The result of this measurement is shown in Fig. 

From this figure one sees that the renormalized Polyakov loop L^en is consistent with a constant in the region 
of lattice spacings investigated. In the same region the unrenormalized Polyakov loop L changes by a factor of 
two. From these facts we conclude that there is no additional ultraviolet divergence in Polyakov loop due to 
non-zero chiral chemical potential. 
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